% This is a script for reading in the phantom files, as an example
Nx = 616; Ny = 485; Nz = 719;

% Load in optical and acoustic phantoms
fid = fopen('./MergedPhantom.DAT', 'r');
phan = fread(fid, 'uint8=>uint8'); phan = reshape(phan, [Nx, Ny, Nz]);
fclose(fid);



